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Q , Abstract. We monitored three transits of the giant gas planet around the nearby K dwarf HD 189733 with the ACS camera 

' on the Hubble Space Telescope. The resulting very-high accuracy lightcurve (signal-to-noise ratio near 15000 on individual 

' measurements, 35000 on 10-minute averages) allows a direct geometric measurement of the orbital inclination, radius ratio 

'. and scale of the system: = 85.68 ± 0.04, Rpi/R, = 0.1572 + 0.0004, a/R, = 8.92 ± 0.09. We derive improved values 

for the stellar and planetary radius, R, = 0.755 ± 0.011 R©, Rpi = 1.154 ± 0.017 Rj, and the transit ephemerides, T,r = 
2453931.12048 ± 0.00002 + n ■ 2.218581 ± 0.000002. The HST data also reveal clear evidence of the planet occulting spots 
. on the surface of the star. At least one large spot complex (> 80000 km) is required to explain the observed flux residuals and 

' their colour evolution. This feature is compatible in amplitude and phase with the variability observed simultaneously from the 

, ground. No evidence for satellites or rings around HD 189733^ is seen in the HST lightcurve. This allows us to exlude with 

■ a high probability the presence of Earth-sized moons and Saturn- type debris rings around this planet. The timing of the three 
transits sampled is stable to the level of a few seconds, excluding a massive second planet in outer 2:1 resonance. 

' Key words, planetary systems ~ stars: individual: HD 189733 - techniques: photometric 
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1. Introduction to first order the orbital inclination and the stellar radius can 

compensate almost perfectly. With ground-based observations. 



The star HD 189733 is transited by a close planetary compan- 



it is extremely difficult to lift completely the degeneracy be- 



?H ion (Bouchy et al. 2005). The brightness of the host star (y=7. 7 , „ ,• , 

rrt ' ; , , • , ,^ ^ , N , , , , tween orbital inclination and stellar radius, because ot system- 

. . , mag), the short period (2.2 days) and the large planet-to-star r- . . i t,, , 

,r ■ ^\ , . , ,- atics induced by the Earth s atmosphere. The spectacular preci- 

radius ratio {Rr,i R, ^0.15) make it the most favourable sys- . i j i i j • • • n j 

A ■^ A A- ^ u f 11 • u A- ^^'^^ reached by space-based transit time series was illustrated 



tem known for detailed studies. In the year following the dis 
covery of its transiting planet, HD 189733 has already been 



with the HST/STIS lightcurve of HD 209458 by Brown et al. 



subjected to many follow-up observations, including extensive ^^^^\u-\u ■ .u in-4 



With an accuracy in the 10 - 10 range, second-order 
features in the transit lightcurve can be detected, such as the 
transit of a satellite of the planet, or the presence of Saturn-type 
debris rings. 



ground-based transit photometry (Bakos et al. 2006, Winn et al. 
2007), measurement of the spin-orbit angle using the Rossiter- 
McLaughlin effect (Winn et al. 2006), measurement and map- 
ping of the infrared flux distribution of the planet with the 
secondary eclipse (Deming et al. 2006, Grillmair et al. 2007, 

Knutson et al. 2007). 2. Observations 

A precise knowledge of the radius of the planet is essen- „, , , . ., , , , . 

. , , ,. ^. , . ,. .f^ We have observed HD 189733 with the Hubble Space 

tial to all applications. From the transit lightcurve, one can ^ 

^, ^ „ ^, ^ , ^, I -, 1 • 1- Telescope during 3 visits of 5 orbits each, using the Advanced 

measure the geometry of the system, hence the orbital inch- ±- o ' o 

^, . r . J • r 1 . 1 1 Camera for Surveys (ACS) in HRC mode with the grism 
nation, the size of the star and the size of the planet. However, j \ j a 
G800L (program GO-10923). In this configuration, the fix?.t- 

Send ojfprint requests to: frederic.pont@obs.unige.ch order spectrum ranges from5500 A to 10500 A with ~ 40 A per 
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pixel. The first visit occured on May 22, 2006 (JD 245877.718 
to 878.092), the second on May 26 (882.115 - 882.408) and 
the third on July 14 (930.946 - 931.240). 

To achieve the highest possible signal-to-noise ratio (S/N), 
we observed with relatively long exposures (25s), allowing the 
ACS detector to saturate by a factor about 5 at the wavelengths 
near 700 nm where it is most sensitive. Prior experience do- 
ing time-series photometry with the ACS showed that with 
proper analysis, Poisson-limited photometry can be obtained 
even with images that are much more saturated than this. This 
is done by integrating over enough pixels to capture essentially 
all electrons generated by the stellar photons, even though they 
may have diffused away from the pixel where they were cre- 
ated. For crowded field photometry, this approach can be prob- 
lematic because of leakage from adjacent star images, but for 
spectroscopy on the bright HD 189733 this is not a problem. 

We obtained extremely high S/N spectrophotometry in a 
time-series mode with a cadence of about 60s. Since the transit 
duration is more than the length of an HST orbit and the target 
star can be seen for less than half of each orbit, three visits 
were required to obtain full phase coverage of the transit. We 
accumulated a total of 675 exposures for the three visits. Each 
visit contains two full HST orbits before the transit, two that 
are all or partly inside the transit, and one that falls entirely 
after the transit. Observations on both sides of the transits are 
necessary for each visit so that clear out-of-transit baselines can 
be established. Also, the first HST orbit often exhibits unique 
systematics as the telescope and instrument stabilize at the new 
pointing. The first orbit of each visit is not utilized in the final 
time series analyses. 

3. Reductions 

Our analysis started with the images provided by the STScI 
pipeline. These are corrected with bias subtraction, dark cur- 
rent subtraction after scaling to individual exposure times, and 
converted from digital units to electrons scaling by the detector 
gain. Retrievals from the STSci archive were made only after 
final biases and darks were available. 

These images are not flat-fielded for grism data. For the 
very similar GO/DD-10441 observations of TrES-1 (Brown et 
al., in preparation) in which HRC grism spectra at high S/N 
were obtained at nearly the same place on the detector, but 
without the complication of saturation and bleeding, it was 
found that use of a wavelength-dependent flat field provided 
no benefit. Facing the extra complication here of saturated data 
and needing to take into account the flat fields as a function 
only of electron origination, not current location, we simply 
have not applied flat fields for these data. For our purposes of 
deriving differential time series over time, in which all of the 
spectra have been taken at nominally the same pointing, and for 
which the extractions are averaged over a large domain (over 
5000 pixels for the global first order spectrum extraction), we 
would not expect much sensitivity to flat fielding errors. If there 
were no motion of the spectra relative to the detector as a re- 
sult of HST guiding errors, there should be no need to flat field. 
These guiding errors, while not zero, are quite small (~ 1/10 



pixel), and at a later step we decorrelate with vectors that re- 
move any small noise resulting from lack of flat fielding. 

3.1. Extraction 

The images consist of a direct image, a first-order spectrum, 
and a partial second-order spectrum, inclined diagonally across 
the detector. The signal of saturated pixels is spread along 
columns (i.e. in the y direction). Analyses start with measure- 
ment of x,y position for the single direct image taken at the 
start of each visit. The calibration given in ACS-ISR 03-07 is 
used to provide a wavelength for each position along the spec- 
trum trace. The global first-order spectrum extraction involves 
a sum over 173 columns that span wavelengths of about 5370 
to 10690 A. The extractions involve sums over columns, rather 
than attempting extractions in the pure cross-dispersion direc- 
tion. For non-saturated data this results in a moderate loss of 
spectral resolution, for saturated data the resolution loss is in- 
herent in the data anyhow. We define an arbitrary data value at 
which to extend the summations along columns, stopping the 
summation the first time the data drop below this value. A num- 
ber of trials are run adopting a range of these arbitrary sum lim- 
iting values, and the one providing highest S/N (after applying 
all the corrections discussed later) is adopted. For HD 189733 
a stopping value of 6,000 electrons works best (peak intensities 
would have been over 10^ in the absence of saturation, thus the 
sums extend to levels of less than 1% relative intensity.) A first 
extraction is done using the above procedure, and intensities 
are normalized by the mean for all points (except first) in orbits 
two and five - this provides an initial light curve. We perform 
a cosmic ray elimination step by sigma clipping data values 
outside the saturated data region (with a buffer of two pixels), 
over the full stack of images in each visit - after normalizing 
out the intensity variations during transit via use of the prelimi- 
nary light curve. After replacing data values affected by cosmic 
rays the normalization is reset. The sigma clipping is done after 
providing a model for the intensity variations based on a linear 
regression with the six external variables to be described be- 
low in the decorrelation section. Cosmic ray elimination is not 
very important for this project that used short exposure times 
and accumulated over 5 x 10^ source electrons each exposure. 
The latter procedure could only catch cosmic rays outside the 
saturated pixel domain. To search for cosmic rays within the 
saturated pixel region we formed stacks of column-summed 
intensities, and performed the same data modelling and sigma 
clipping step as above. This resulted in trivial improvements for 
2 of 3 visits, and a minor loss for the third. We have kept this 
step although it neither helped nor harmed the overall result to 
any significant extent. 

For a specific application (see Section 14.5b . we also ex- 
tracted sub-sections of the spectrum in 500 A wavelength in- 
tervals, restricting the extraction to specific ranges in y. 

3.2. Decorrelation 

It is important to establish by how much the spectrum moves 
relative to the detector exposure-to-exposure, and at the same 
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time establish records of other characteristics such as measures 
of the cross-dispersion width and the rotation of the spectrum 
relative to the detector. Such external vectors become candi- 
dates for performing decorrelations with the intensity vectors. 

For each time step a one-dimensional Gaussian was fit to 
all columns in the first- and second-order spectra for which (a) 
the core does not saturate, (b) the column is not next to one that 
saturates, and (c) the column intensity sum reaches 10^ elec- 
trons. The results of these fits are then averaged together. After 
subtracting the mean this provides time series records of: (1) 
mean y-position of the spectra, (2) a measure of the spectrum 
width - Gaussian sigma, and (3) a rotation formed by differenc- 
ing mean y-position in first and last 10 columns, which for these 
data are separated by about 65 columns. A measure of the spec- 
trum shifts in x are derived by forming a one-dimensional vec- 
tor of column-summed intensities over x, averaged over all the 
exposures in a visit. Pixel-by-pixel in x the intensity derivative 
in X is formed, and note is taken of the region (columns 2-15 
in practice) that has both decent signal and strong derivatives. 
Then averaged over these 14 columns the relative intensities 
in individual spectra in combination with the dl/dx derivative 
provide delta-jc shifts for each exposure. The spectrum shifts 
in X with a characteristic pattern each orbit by about 0. 1 pixel 
peak to peak. A somewhat larger offset occurs in y, with both a 
long-term drift and within orbit variations present. The cross- 
dispersion width becomes more variable with each successive 
orbit. Rotation maintains a characteristic pattern each orbit. 

Decorrelations are performed by evaluating a multi-linear 
regression on these vectors simultaneously: x-position, y- 
position, width, rotation, HST orbit phase, and overall time. 
No attempt is made to account for the fact that several of these 
external vectors show mutual correlations. The regression is 
based on a least-squares fit using all of the points outside of 
transit, ignoring orbit 1 . Then the full set of data points is di- 
vided by the fit to provide corrections for noise that correlate 
with these external vectors. 

The decorrelations generally remove noise with character- 
istic time scales smaller than the HST orbit, plus an overall 
linear trend. For the average number of 130 data points used 
to establish the decorrelation fits in each visit, the time series 
standard deviation improves from a geometric mean over the 
three visits of 277 • 10"^ before, to 67 • lO"*" after decorrelation. 
The mean absolute value of individual linear correlation coef- 
ficients are: x-position (0.22), y-position (0.48), width (0.35), 
rotation (0.53), orbit phase (0.32), and overall time (0.84). 

No additional sigma-clipping of the data has been per- 
formed. Exposure 65 of visit 02 is affected by a strong cosmic 
ray and was not used. 

3.3. Systematics after decorrelation 

Since a primary objective of this program is to measure devia- 
tions from the transit shape at the level of a fraction of a milli- 
magnitude, it is crucial to estimate the level of systematics that 
could remain in the data after decorrelation. One method is to 
use decorrelation trials in which some out-of-transit orbits are 
not used in the regression to estabhsh the decorrelation. Here 



there are two visits, the first (VI) and second (V2) that have 
two orbits either after (VI) or before (V2) the transit. The ex- 
ercise is to step through not using either orbits 2 or 3 of V2, or 
orbits 4 or 5 of VI. The residuals are then compared to those 
in which all of the data were used. The results were very clear 
In the four orbits (total) not used in decorrelations the largest 
deviations remained under 10 with means typically an order 
of magnitude smaller still. In no case did any residual show de- 
viations as large as those seen multiple times for orbits inside 
the transits. This implies that the variations are real and not due 
to instrumental systematics. 

A second method, this time to check if some real signal 
could have been removed by the decorrelation, is to insert fea- 
tures in the data before decorrelation, and check how they are 
recovered. Decorrelation trials in which signals similar to those 
in V2 (rise from 0.0 to +Q.5 mmag in about 10 points followed 
by a drop back to 0.0 in 10 points) are added to the time series 
before decorrelation. We chose to add these to either orbits 2 or 
3 of V2, or 4 or 5 of VI (one at a time). In 4 of 6 total cases this 
triangle-shaped perturbation was added to centers of orbits, in 
the other 2 it was added near the end of orbits. In each case 
the signal is preserved through the decorrelation step at the 80- 
90% level. It is not surprising that there is a slight reduction in 
height resulting from having added in a positive signal. 

The above two tests show that: (a) the out-of-transit orbits 
stay less noisy, and show no systematics even when not used to 
establish the decorrelation, and (b) if real systematic deviations 
had been present that last only a fraction (~l/3rd) of an orbit, 
these would only be affected at the 10-20% level by decorrela- 
tion. This adds significantly to believing that the residual fea- 
tures seen during the in-transit HST orbits are real, and likely 
reflect flux inhomogeneities on the surface of the star. 

We also tested that our results are not sensitive to the size of 
the extraction box used, by repeating the analysis with a limit 
at 2,000 electrons (instead of 6,000). 

3.4. Absolute flux 

Previous experience with ACS and HST show that the variabil- 
ity of ACS over a few weeks is about 0.003 to 0.004 in frac- 
tional flux calibration. Much of the deviation is correlated with 
known variation of focus. After correcting for these variations, 
an epoch-to-epoch precision of ~< 0.002 can be obtained. 

The visit-to- visit changes of total intensity are, based on the 
mean of orbits 2 and 5 without the first point in orbits: -1-0.005 
for visit V2 relative to visit VI, and -hO.007 for V3 relative to 
VI. The Poisson limit on these would be less than 5- 10"^. There 
is a caveat that needs to apply to V3 - this is after the switch to 
Side 2 electronics in early July 2006. There could be a minor 
change in operating temperature. We don't expect a change, but 
there is a potential zero-point difference for visit V3. 

4. Results 

4.1. Lightcurve 

The lightcurve decorrelated from dependence on instrumental 
parameters is given in Table[T] and shown in Figure[T] The data 



4 



F. Pont et al.: Hubble Space Telescope times-series photometry of HD189733 



1 ^ 




Phase 



Fig. 1. DecoiTelated lightcurve phased to P=2.2 18583 days, with best-fit model transit curve. Top: flux after external parameter 
decoiTelation as a function of phase; Bottom: residuals around the best-fit transit model. Light blue for the first visit, dark blue 
for the second and green for the third. Open symbols indicate data affected by Features A and B (see text), not used in the fit. 
Dotted lines show the +10"^ level. 



Date 


Flux 


2453877.71866 


0.999421 


2453877.71935 


0.999318 


2453877.72005 


0.999580 


2453877.72074 


0.999518 


2453877.72143 


0.999498 


2453877.72213 


0.999603 


2453877.72282 


0.999592 


2453877.72352 


0.999637 


2453877.72421 


0.999532 


2453877.72491 


0.999610 



Table 1. Photometric time series (electronic table) 



in the first orbit of each visit, not used in the analysis, are not 
shown. 

The standard deviation of the individual data points outside 
the transit is 67 x 10"*. Expressed in terms of signal-to-noise, 
this corresponds to 15000:1. The photon noise is 45 x 10 
The standard deviation on 10-points average is 28 x 10"*, cor- 
responding to a S/N of 35000: 1 . 



No out-of-transit structure is visible in our lightcurve above 
the level of the residual instrumental systematics. 

Transit lightcurves are expected to be perfectly symmetri- 
cal to better than the 10""* level. Two asymmetric features are 
apparent in our data during the transit, one at the end of the 
second orbit of the first visit, the other at the beginning of the 
third orbit of the second visit. We label the first "Feature A" in 
our analysis, and the second "Feature B". 

Feature A is much too large to be explained by an instru- 
mental effect, and no special behaviour of the telescope was 
observed at the corresponding time. The flux increase during 
Feature A is also seen in the zeroth-order image on the CCD. It 
is accompanied by a detectable change of spectral distribution, 
which suggests an explanation in terms of the transiting planet 
occulting a cool spot on the surface of the star (see Section l43] l. 

Feature B is also larger than instrumental effects, and does 
not correlate with any of our external instrumental parameters. 
Its less regular shape and the fact that any colour effect is be- 
low the noise level indicate that, in principle, an explanation in 
terms of instrumental noise cannot be entirely excluded. Based 
on our experience with previous HST high-accuracy times se- 



F. Pont et al.: Hubble Space Telescope times-series photometry of HD189733 



5 



ries, and the simulations of Section 13.31 we believe, however, 
that Feature B is also real. 

We do not use Features A and B in the analysis of the 
lightcurve in terms of planetary transit. They are treated sep- 
arately in Section l43] 

4.2. Stellar activity and variability 

HD 189733 is an active star, variable to the percent level. It 
is listed in the Variable Star Catalogue as V452 Vul. A chro- 
mospheric activity index of S - 0.525 has been measured by 
Wright et al. (2004). Activity -related X-ray emissivity has been 
measured by both EXOSAT and ROSAT, activity-related ra- 
dial velocity residuals of 15 m ' were reported by Bouchy et 
al. (2005). Winn et al. (2007) have measured the photometric 
variability of this star extensively and confirm variability at the 
percent level, compatible with an explanation in terms of tran- 
sient spots modulated by a rotation period of Prot ~ 13.4 days. 
Moutou et al. (2007) measure strong activity in the Call line 
and infer a strong magnetic field with a complex topology from 
spectropolarimetric monitoring. The explanation of Features A 
and B in terms of starspots is therefore natural. The presence 
of large starspots is also confirmed by an observing campaign 
on this object by the MOST satellite (CroU et al., in prep.). The 
MOST data yield an improved rotation period of 11.8 days. 

The Winn et al. (2007) photometry is contemporaneous 
with our HST data. Our absolute measurements are placed 
within the context of the ground-based monitoring in Fig. |2] 
with an arbitrary zero-point shift. The HST data is in agreement 
with the periodic variation seen in the long-term lightcurve. If 
we interpret this variability in terms of starspots moving in and 
out of view with the rotation of the star, then the third visit 
occurs near the brightest point - with less star spots visible - 
and the first visit with a 0.007 dimming due to starspots. The 
phasing and amplitude of features A and B are perfectly com- 
patible with an explanation in terms of the planet occulting part 
of the starspots responsible for the photometric variation (see 
Sectionl43]i. 

Before fitting a transit signal, we correct for the variations 
of the total stellar luminosity due to the presence of starspots. 
Outside of features A and B, the planet crosses a spot-free re- 
gion, therefore a region slightly brighter than the average over 
the stellar disc, which includes the spots. This is a tiny correc- 
tion of the scaling between transit depth and radius ratio (of the 
order of 2 x lO""* in flux). Nevertheless, to the level of the ac- 
curacy of the HST lightcurve, it makes a significant difference 
and must be accounted for. We use the absolute flux differences 
measured in Section[3l4l 



4.3. Transit signal 

A transit light curve computed with the Mandel & Agol (2002) 
algorithm was fitted to the light curve, with a downhill simplex 
algorithm (Press et al. 1992). Features A and B were removed 
with cuts from JD=877.875 to the end of orbit 3 of visit 1, and 
from the beginning of orbit 3 of visit 2 to JD=882.333. The 
limb-darkening coefficients were left as free parameters, us- 
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Fig. 2. Photometry of HD 189733 from Winn et al. (2007; 
dots), with the timing and mean flux level of the three HST 
visit indicated (shifted by an arbitrary constant; stars). 



Measured in this paper 
Period [days] 
Transit epoch [JD] 
Radius ratio Rpi/R, 
System scale a/R, 
Impact parameter a cos i/R, 
Orbital inclination ; ["] 
Host stai M"^/R ratio 
Visit 1 transit epoch [JD] 
Visit 2 transit epoch [JD] 
Visit 3 transit epoch [JD] 

Inferred 

Stellar radius [Rq] 
Planetary radius [R j] 



2.218581 ±0.000002 
2453931.12048 ±0.00002 
0.1572 ±0.0004 
8.92 ±0.09 
0.671 ±0.008 
85.68 ±0.04 
1.242 ±0.012 
2453877.87448 ±0.00005 
2453882.31171 ±0.00005 
2453931.12048 ±0.00003 



0.755 : 

1.154: 



lO.Oll 

:0.017 



Adopted (from Bouchy et al. 2005) 

Stellar mass [Mq] 0.825 ±0.025 

Planetary mass [Mj] 1.15 ±0.04 

Table 2. Parameters of the transiting system HD 189733. 



ing the four-coefficients non-linear expression of Claret (2000). 
They converge to a i_4 = (0.506,-0.727, 1.345, -0.519). These 
values are compatible with the ones for Tgff = 5000K, 
[Fe/H]=0.0 and logg = 4.5 for the R filter. 

The resulting best-fit transit parameters are shown in 
Table |2] The best-fit theoretical transit curve and the residuals 
are shown in Fig.[T] 

The standard deviation of the residuals around the model 
lightcurve in the transit are 129x10"^. This is larger than the 
photon noise, and is due to the presence of systematic fluctua- 
tions on the scale of 10 even outside Features A and B. These 
fluctuations may be due to systematics not perfectly corrected 
by the decorrelation, or to structures on the surface of the star 
smaller in scale than those responsible for Features A and B. 

Including parts or all of Features A and B in the fit, or 
changing the limb darkening coefficients within reasonable 
limits, does not modify the resulting system parameters beyond 
the 1 % level, and therefore has no influence on the astrophysi- 
cal applications. 

The photon-noise uncertainties on the transit parameters 
are extremely small, and the actual uncertainties are dominated 
by correlated flux residuals. The amplitude of the residuals and 
the tests in Section [3] show that these eff'ects are of the order 
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of 10""* at most over the relevant timescales. We therefore esti- 
mate the uncertainties on the transit parameters by allowing the 
data to move, as a whole, by a maximum of 10"'* in flux from 
the model, following the approach of Pont, Zucker and Queloz 
(2006) to derive error intervals in the presence of correlated 
noise. These values can be considered upper limits for the un- 
certainties, since the systematics are unlikely to be exactly in 
phase with the diff'erences caused by changing the transit pa- 
rameters. 

The uncertainties due to these fluctuations, in turn, are 
dominated in the derivation of the physical parameters of the 
system by the uncertainty on the stellar mass and radius, due 
to the RM^^^^ degeneracy of transit systems. As an order-of- 
magnitude indication, the photon noise uncertainty on the ra- 
dius ratio is ~ 0.01%, the uncertainty due to systematics is ~ 
0.1 %, and the uncertainty due to RM^^^^ degeneracy is ~ 1% 
(See Sectionl43]i. 

4.4. Transit timings 

Table |2] also gives the values of the central epoch of the three 
transits covered by our HST data, when fitted independently 
with the other parameters fixed at the values given by the com- 
bined data. 

We compared these values with previous ground-based 
mesurements (see Figure[3j. Deviations of transit timings from 
strict periodicity can reveal the existence of other planets in 
the system (e.g. Agol et al. 2005, Holman & Murray 2005). 
However, the relation between transit timing variations and 
the orbital elements of the unseen perturber is very complex, 
and non-detections exclude only restricted regions of param- 
eter space. Periodicity at the level of a few seconds excludes, 
for instance, the presence of an Earth-mass body in 2:1 reso- 
nance with the giant planet (see figure 5 of Agol et al. 2005). 
The HST transit timings are compatible with previous deter- 
minations, while being much more precise. The periodicity of 
the transits is stable at the level of a few seconds over the three 
HST transits (the residuals relative to the periodic solution are 
5+4 s, -1 +4 s and 0±3 s for the three visits respectively). This 
strongly suggests that the excess seen in the transit timing 
for ground-based data is due to unrecognised systematics, as 
suggested by Bakos et al. (2006), rather than real orbital varia- 
tions. 

The only transit timing of precision comparable to the HST 
is that of Knutson et al. (2007), the last point on Figure [3] It 
is incompatible with the HST timings at the ~ 4cr level. This 
could be a first indication of transit timing variations, of the 
order of half a minute. Given the difficulty of estimating sys- 
tematic errors at this level of accuracy, more data is needed to 
confirm this tentative indication. An intense monitoring of HD 
189733 with the Spitzer space telescope is currently under way. 

4.5. Moons, rings, and spots 

No evidence is seen in the flux residuals for the characteris- 
tic signature of: (a) a transiting satellite or second planet (b) 
planetary rings (c) planet or star oblateness (d) transit timing 
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Transit timing [JD-2450000] 

Fig. 3. Transit timing residuals for HD 189733, from Bakos et 
al. (2006), Winn et al. (2007), Knutson et al. (2007) and this pa- 
per, relative to our best-fit ephemerides. The three HST transits 
are indicated by the larger squares. The dotted lines delienate 
our 1 - cr uncertainty interval. 
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Fig. 4. Residuals around the best-fit transit model. Open sym- 
bols indicate data during the transit, crosses data attributed to 
Features A and B. The lines show the signals expected from an 
additional transiting 0.8 Earth-radius body (solid), from plan- 
etary rings (dotted), and from transit timing variations of 20 
seconds (dashed). 



variations. The residuals around the best-fit transit lightcurve 
are compared to the expected shapes for three effects in Fig.|4] 
(the eff'ect of planet oblateness resembles that of rings). 

The two asymmetrical features observed during the transits 
do not have the correct shape to be explained in terms of the 
transit of a second body in the system. They do not occur at 
the phase expected for such effects as the presence of rings, 
non- sphericity of the star or planet. Feature A and B are best 
explained as due to the passage of the planet in front of cooler 
regions - starspots or spot complexes - on the surface of the 
star itself. 



F. Pont et al.: Hubble Space Telescope times-series photometry of HD189733 



7 



a. 





1 

\ 


1 


1 1 


1.2 


\ 






1 




X 




- 




- 


0.8 


- , 1 , 


, 1 , 


V i-J 1 1 r 




6000 


7000 


8000 9000 



Wavelength [A] 

Fig. 5. Slope of the residual flux rise caused by "Feature A" in 
the time series of 500 A spectral bins, divided by the slope of 
the lightcurve integrated over all wavelengths. The lines show 
the expectation for a spot 1000 K cooler than the surface of the 
star (short dashes), 100 K cooler (long dashes), and for system- 
atic noise without a colour signature (dots). 

Feature A, to first order a linear rise in flux, can be well 
explained by the planet occulting part of a large cooler region 
on the star. Feature B is less regular and has lower S/N. Its 
exact shape is sensitive to uncertainties on the limb darkening 
parameters and the absolute flux difference between the HST 
visits. 

For the duration of Feature A, the planet moves ~ 80 000 
km across the face of the star The orbital and rotation axes 
of the system are almost aligned (Winn et al. 2006) and the 
rotation velocity of the star is long compared to the duration of 
the transit (Winn et al. 2007), so that this distance corresponds 
to a similar projected distance in longitude on the surface of the 
star This gives the minimal extension in longitude of the spot 
responsible for Feature A. 

The minimum latitude extention of the spot is given by the 
total flux increase of Feature A. If the spot is much cooler than 
the rest of the star surface, the flux it emits is negligible, and it 
can be considered a dark feature. In that case, the width needed 
to produce the observed flux increase is 12 000 km. At the other 
extreme, the occulted portion of the spot cannot be larger than 
the diameter of the planet, R~ 165 000 km. To produce the ob- 
served drop in flux, a spot of this size will need to have a tem- 
perature less than 100 K degrees cooler than the rest of the star. 

There is one additional constraint in the data to support the 
spot interpretation. Since a starspot is cooler than the suiTound- 
ing regions on the star's surface, the occultation of a spot will 
introduce a wavelength dependence in the flux rise. The flux 
rise will be steeper in the blue, where the temperature depen- 
dance is sharper, than in the red. The magnitude of this effect 
will depend on the effective temperature of the spot compared 
to the rest of the star 

We computed the flux in our spectra in 500-A wide spec- 
tra bins, from 5500 to 9500 A (decorrelated for systematics 
with the procedure explained in Sect. |3]l. For each bin, we fit- 
ted a transit model lightcurve with the parameters fixed as per 
Table |2] leaving only the limb-darkening coefficients free. We 
then fitted a linear slope to the flux rise caused by Feature A. 



We divided this slope by the slope measured on the integrated 
lightcurve. Figure|5]shows the value of this slope as a function 
of the central wavelength of the spectral bin. We compare these 
to the expectations for a signal caused by occulting a cooler re- 
gion (dashed lines in the figure), using the blackbody approx- 
imation for the spectral distribution of flux in the stellar atmo- 
sphere inside and outside the spot and setting T=5000 K for the 
stellar surface outside the spot. The data is clearly compatible 
with the spot scenario, and a colour-independent systematics 
(dotted line in the figure) can be excluded. 

In principle, this procedure could also allow us to measure 
the temperature of the spot. But in this range of temperatures 
and spectral coverage, a bigger, warmer spot and a smaller, 
cooler spot produce very similar signatures, that are within the 
observational uncertainties. 

In summary. Feature A is compatible with a complex of 
cool spots extending at least 80 000 km in longitude and 12 000 
km is latitude, in a relatively uniform manner, with an overall 
effective temperature much lower than the rest of the star It is 
also compatible with a larger, more circular spot, with a smaller 
temperature difference, down to about AT" <~ 100 K for a spot 
as wide or wider than the planet. 

The same procedure was repeated for Feature B. However, 
in this case, the total amplitude of the effect is too small, and 
any colour signature is dominated by the uncertainties. Also, 
the more complex shape of the signal is compatible with sev- 
eral arrangements of spots. Therefore, while we posit that the 
flux variations corresponding to Feature B are due to the oc- 
cultations of activity-related features on the surface of the star, 
we do not attempt to model this feature in detail. In the limit of 
large temperature difference (dark spots), the size of the spots 
responsible for Feature B would have to be equivalent to a cir- 
cle of 9 000 km radius (more than one spot would be necessary 
to reproduce the observed lightcurve shape). 

Figure |6] illustrates the geometric configuration of the star 
and planet as seen from Earth during the three HST visits, giv- 
ing a schematic indication of the position and size of the spots. 
The path of the planet is shown only during HST observations 
(at other times the HST line of sight is blocked by the Earth). 
A rotation period of 1 1 .8 days is used to rotate the star between 
the different visits. The rotation of the star during the transits 
is very small (about 2 degrees), and can be neglected. There 
is a global north-south degeneracy. The longitudinal extension 
of the spots is well constrained, but not their exact sizes and 
latitudinal position relative to the planet. The stellar rotation is 
aligned with the planetary orbit, as measured by Winn et al. 
(2007) using the Rossiter-McLaughlin effect. 

The spot responsible for Feature A is not encountered again 
during the two other visits. However, the spot complex caus- 
ing Feature B is transited again during egress in visit 3. This 
may explain the behaviour of the residuals at that point in the 
lightcurve (see bottom panel of Fig. Altogether, the planet 
transits about 10% of the star's surface during our observations, 
and occults starspots about half this time, suggesting that spots 
of different sizes are abondant on the stellar surface. 
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Fig. 6. Schematic configuration of tlie position of the star and planet during the tliree HST visits. Longitude and latitude hues are 
drawn on the star at 30 degrees intervals. Longitude zero is defined as the longitude pointing towards the Earth at the epoch of 
the first transit. The arrows indicate the direction of the stellar rotation and planet orbital motion. The path of the planet during 
the HST observations is shown for the three visits. Darker features on the star sketch the configuration for the spots affecting the 
HST lightcurve. 



4.6. Host Star and planet 

The radius of the primary and its mass to the one-third power 
are degenerate in transit lightcurves. Second-order effects that 
would allow breaking this degeneracy, such as the light travel- 
time effect, are far too small to be detectable in the presence of 
any realistic amount of photometric noise or stellar variability. 
Therefore, the M '^^ R degeneracy must be lifted with recourse 
to stellar evolution models. 

Fortunately, HD 189733 is a well-observed, low-mass star, 
and evolution models allow for only a small range of masses 
compatible with its observed visible and infrared magnitude, 
parallax, and spectroscopic temperature. Its Hipparcos paral- 
lax of 51.94 +0.87 mas implies a distance of 19.25 +0.32 par- 
sec, and absolute magnitudes of Mv=6.25 in the visible and 
M/f =4.07 in the infrared. These luminosities can be combined 
with iheM'^'^R- 1 .246 ± 0.0 1 2 (in Solar units) value from the 
HST transit shape and confronted to stellar evolution models. 
Figure |7] shows the error box for HD 189733 in these two pro- 
jections of parameter space compared to the position of Padova 
stellar evolution models (Girardi et al. 2002). The data are co- 
herent with the models and compatible with a mass in the 0.80- 
0.85 Mq range for HD 189733. Similar results are obtained us- 
ing the Baraffe et al. (1998) low-mass star models. We therefore 
confirm the value of M = 0.825 + 0.025 Mq used in Bouchy et 
al. (2005) for the mass of the host star. 

It is interesting to note that low-mass eclipsing binaries 
with well-determined mass and radius are not compatible with 
the model tracks in this range of mass (Ribas et al. 2006), un- 
like HD 189733. This is a strong indication confirming the 
suspicion that close eclipsing binaries strongly influence each 
other, and are therefore not good calibrators for the mass-radius 
relation of single stars. 

We adopt M - 0.825 + 0.025Mo to lift the mass-radius 
degeneracy, which implies Rs,ar - 0.753 + 0.01 IRq, and Rpi = 
1.147 + 0.0177;y (with Rj = 71500 km). These values are listed 
in Table |2] 

The planetary mass remains unchanged at Mpi - 1.15 + 
0.04 from Bouchy et al. (2005). This corresponds to a planetary 
density of 822 kg m"^. 



5. Discussion 

The ACS spectrophotometric time series presented in this ar- 
ticle shows the capacity of HST to obtain lightcurves of ex- 
treme accuracy for bright stars. It echoes the HST lightcurve 
of Brown et al. (2001) acquired with the STIS spectrograph 
for HD 209458. The S/N for HD 189733 is even higher be- 
cause of the larger wavelength domain covered (at lower spec- 
tral resolution), the fact that HD 189733 is slightly brighter in 
the wavelength interval sampled, and the higher throughput of 
ACS/HRC. 

The single-point standard deviation outside the transit is 
67 X 10 '', for a sampling of one point per minute. Correlated 
systematics left after decorrelation of instrumental parame- 
ters such as tracking drift and focus change have an ampli- 
tude well below the 10 level. The standard deviation of the 
mean of measurements over 10 minutes is 28 x 10"^. By com- 
parison, systematics in high-precision ground-based photom- 
etry are typically of the order of one part per thousand, or 
1000 X 10"^ 

These data allow a very precise measurement of the plane- 
tary radius to be obtained. This is also due to the fact that the 
constraints on the mass of the star are good, because it is a low- 
mass star with a known parallax from Hipparcos. 

The HST/ACS lightcurve permits to lift the degeneracy be- 
tween primary radius and orbital inclination in the transit sig- 
nal. Bakos et al. (2006) show how the data of a ground-based 
photometric campaign, even with several telescopes and in sev- 
eral filters, can lead to misleading conclusions in this respect. 
The Bakos et al. data apparently indicated a much smaller pri- 
mary (Rt = 0.68 + 0.02 Rq), but the authors concluded that this 
was due to unrecognized systematics in the ground-based pho- 
tometry. We reinforce this conclusion. Our results, on the other 
hand, confirm the values of Winn et al.(2007), with smaller 
uncertainties. We conclude that Winn et al. shows the proper 
amount of data necessary to lift the radius-inclination degener- 
acy from the ground with confidence, for such a deep transit: 8 
complete transits covered in excellent conditions. This gives a 
yardstick against which to evaluate similar ground-based cam- 
paigns. 
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Fig. 7. Constraints on HD 189733 from the absolute visible and infrared magnitude, and the M'^^ 7? ' ratio, compared to Padova 
stellar evolution model isochrones. 



More recently, Baines et al. (2007) have measured the ra- 
dius of HD 189733 direcdy with the CHARA interferometer, 
obtaining /?» - 0.779 + 0.052/?o. This value has higher un- 
certainties but is completely independent of any assumption of 
mass, and is therefore a useful check to the determinations from 
the geometry of the transit. It is compatible with our value. 

From the 8-;U transit lightcurve measured with the Spitzer 
space telescope, Knutson et al. (2007) obtain Rpi/Rsm = 
0.1545 + 0.0002. This is significantly lower than our value. 
Starspots not crossed by the planet during the transits can con- 
tribute to this difference. These spots will lower the surface 
brightness of the part of the star not occulted by the planet, 
therefore increasing the flux drop and the apparent radius ratio. 
This effect is smaller in the infrared, because the spot contrast 
gets smaller. Another possible explanation is that the effective 
transit radius of the planet is indeed larger in the visible and 
near infrared than in the mid-infrared, due to the transmission 
spectrum of the planetary atmosphere. These topics are beyond 
the scope of this paper. 

We here quantify our non-detection of moons or rings in 
terms of upper limits to any satellite or ring system around HD 
189733b. 

We have an upper limit of 10 to the depth of any unrecog- 
nised second eclipse in our data. This corresponds, given the 
size of the primary star, to a body of 5000 km (0.8 Earth radii). 
Unless the satellite spends each of the three transits that we 
have sampled either in front or behind the planet - an extremely 
unlikely configuration - we conclude that a moon larger than 
0.8 Rg) around HD 189733b is excluded by the data. Figure E] 
shows the signal such a moon would cause, superimposed to 
the residuals during the three HST visits. 

We also exclude a ring system that would cause a flux drop 
larger than 10 around the time of transit. Since the planet is 
orbiting very close to the star and is subject to strong tides, we 
assume that any ring system would be aligned with the orbital 



plane (/ - 85.72). In that case, a Saturn-type ring system at 
about 1 Rpi from the star's surface will need to be less that 
120 r km wide, where t is its optical depth (if Saturn's ring are 
any guide, a ring system should be thin enough that the ring's 
thickness be negligible compared to its projected width even 
for the low inclination of HD 189733, so that this result is vaUd 
for any viewing angle of the rings except exactly edge-on). For 
comparison, the Saturn main ring system is 40 000 km wide 
with a high optical depth (~ 0.3), and the Jupiter ring system 
is also thousands of km wide but has a very small optical depth 
(10"^ to 10"^). We can therefore exclude a major debris ring 
around HD 189733b (discounting the unlikely coincidence of 
an exact alignement with the line-of-sight). 

We note that altogether these limits would not have allowed 
us to detect any of Jupiter's moons, and Jupiter's weak ring sys- 
tem. Therefore, the absence of moon and rings in HD 189733 
shown by our data must be understood as showing the absence 
of major. Earth-size moons and Saturn-type debris rings. 

The fact that spot-like features afi'ected a large part of the 
HST transits suggests that activity-related features are abun- 
dant on the surface of HD189733. The fliree HST visits to- 
gether occulted slightly less than 10% of the star's surface, re- 
vealing one spot at the 10"^ level and one at the 4 ■ 10"^ level. 
This is in good order-of-magnitude agreement with the vari- 
ability observed by Winn et al. (2007), which changes from 
negligible to 1% over a few rotation cycles. 

A coherent picture of the active surface of HD189733 
emerges from the combination of the small-scale information 
from HST data and the large-scale variability monitored by 
Winn et al. (2007). The timing of the sudden flux rise seen at 
the end of the first HST visit ("Feature A") is compatible with 
the planet occulting part of a large spot responsible for a sig- 
nificant fraction of the modulation of the lightcurve. Between 
the first and the second visit, the luminosity of the star has in- 
creased by 0.5%, and about a third of the stellar rotation period 
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has elapsed, giving the spot (or spots) enough time to rotate 
from sUghtly east of the center of the star to the hidden side. 

In conclusion, our results demonstrate again the potential 
of HST to gather extremely precise lightcurves for bright stars. 
We derive improved new measurements of the characteristics 
of this important system, that provide very tight constrains 
for the models. Our lightcurve did not show any evidence for 
Earth-like satellites or rings around the planet HD 189733, nor 
transit timing variations indicative of an unseen second planet. 
The lightcurve during the transit was strongly affected by the 
spots and active regions crossed by the planet, providing a rare 
glimpse of the small-scale geometry of spots on a star other 
than the Sun. 
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